Transistor-Like Behavior of a Bose-Einstein Condensate in a Triple Well Potential. 
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In the last several years considerable efforts have been devoted to developing Bose-Einstein Con- 
densate (BEC)-based devices for applications such as fundamental research, precision measurements 
and integrated atom optics. Such devices capable of complex functionality can be designed from 
simpler building blocks as is done in microelectronics. One of the most important components of 
microelectronics is a transistor. We demonstrate that Bose-Einstein condensate in a three well po- 
tential structure where the tunneling of atoms between two wells is controlled by the population in 
the third, shows behavior similar to that of an electronic field effect transistor. Namely, it exhibits 
switching and both absolute and differential gain. The role of quantum fluctuations is analyzed, 
estimates of switching time and parameters for the potential are presented. 
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I. INTRODUCTION 

Recent experimental realizations of atom optical de- 
vices such as atomic waveguides, beamsplitters 01 [USE! 
on-chip BEC sources and conveyor belts 0, |f| opens a 
way for development of more complex devices such as, 
e.g., BEC-based interferometers 0, 0]. On-chip inte- 
grated cold atom circuits capable of complex function- 
ality can be constructed from simpler building blocks as 
it is done in microelectronics to find applications in fun- 
damental physics, precision measurements and quantum 
information technology. 

One of the most important components of a microelec- 
tronic circuit is a transistor. In this paper we present 
a BEC-based device which will be subsequently called a 
BEC transistor or an atom transistor. It enables one to 
control a large number of atoms with a smaller number 
of atoms and demonstrates switching and both differen- 
tial and absolute gain thus showing behavior similar to 
that of an electronic transistor. The device is not opti- 
mized for performance but is arguably the simplest pos- 
sible geometry showing behavior reminiscent of a tran- 
sistor. This makes its experimental realization relatively 
easy with existing atom chip techniques. 

The BEC transistor uses a Bose Einstein condensate 
in a triple well potential, as shown schematically in Fig 
^ hi fact, Fig. ^ refers to two subtly different possible 
experimental realization of the device. In the trapped 
configuration, the BEC is confined in all three dimen- 
sions in the potential wells. The wells are allowed to 
interact for time interval T. This is done either spatially 
bringing them together and separating apart after time T 
of changing the shapes of the potential wells so that the 
interaction is suppressed after time T. In the waveguide 
configuration, the potential wells of Fig.^represent three 
guides that converge, run parallel to each other for dis- 
tance L and then diverge. The interaction time T = L/v 



in this geometry is determined by the speed of flow v of 
the BEC in the guides. In the following for definitcness 
we will use the terminology appropriate for the trapped 
configuration. 

The BEC transistor is similar to an electronic field ef- 
fect transistor. The left well behaves like the source, the 
middle as the gate, and the right well is equivalent to the 
drain. If there are no atoms in the middle well, practi- 
cally no atoms tunnel from the left into the right well, 
as shown in Fig. [IJi. A small number of atoms placed 
into the middle well switches the device resulting into 
the strong flux of atoms from the left well (the source) 
through the middle and into the right well as shown in 
Fig. \T]p. Increasing number of atoms in the middle well 
increases the number of atoms that tunnel into the right 
well. Parameters of the triple well structure are chosen so 
that the number of atoms having tunneled into the right 
well at the end of the interaction period is much larger 
that the number of atoms in the middle well. In the sub- 
sequent sections we will show that the BEC transistor 
exhibits both absolute and differential gain. 

The physics of operation of the BEC transistor is based 
on atom-atom interactions and appropriate design of the 
potentials. The chemical potential of the left well is cho- 
sen to be nearly equal to the ground state energy level 
of the empty right well (in Fig. ^ we make them equal). 
The ground state energy of the empty middle well is cho- 
sen to be considerably lower than that in both the left 
and the right wells. Placing atoms in a well raises the 
value of chemical potential due to atom-atom interac- 
tions. Parameters of the potential wells are chosen so 
that the chemical potential in the middle well is con- 
siderably more sensitive to the change in the number of 
atoms in the well than is the case for the left and right 
wells. 

When the middle and right wells are initially unpop- 
ulated, tunneling of atoms from the left to the middle 



well is blocked because of the energy mismatch as shown 
in Fig. If some amount of atoms is placed into the 
middle well, the atom-atom interactions will increase the 
energy of the atoms in the middle well. When the chem- 
ical potential in the middle well becomes nearly equal to 
that in the left and right wells, the device switches and 
atoms become able to tunnel from the left through the 
middle into the right well as shown in Fig. ^1. 

Using atom-atom interactions to block tunneling in a 
double- well structure is often referred to as self trapping. 
This effect was first described in Ref. |9| . If a condensate 
is placed in one of the two weakly-coupled spatially sep- 
arated potential wells with matched energy levels, it can 
oscillate between the wells by linear quantum tunneling. 
However, due to atom-atom interactions, the tunneling 
is blocked when the number of atoms in the condensate 
exceeds some critical value. This suppression is due to 
the fact that interactions increase chemical potential of 
the atoms in the occupied wells and introduce nonlin- 
ear energy mismatch. Self trapping has been analyzed 
for a large number of systems including asymmetric dou- 
ble well potentials |10| and symmetric three well systems 
PH. It has also been observed experimentally for atoms 
in a one dimensional optical lattice 01 ■ 

The quantum state of two trapped Bose-Einstein con- 
densates in a double well potential has been analyzed in 
Ref. It has been shown that when the two wells 

are separated and the interaction between the atoms is 
repulsive, the lowest energy state is fragmented, which 
means that the coherence between the atoms in each well 
is lost. The dependence of this fragmentation on the 
splitting rate and physical parameters of the potential 
has been analyzed in Refs. fl4| and 0|. The visibility 
of interference fringes after splitting of a condensate with 
both attractive and repulsive interactions was analyzed 
in Ref. [l6|. who showed a decrease in quantum noise in 
the case of attractive interactions. The quantum dynam- 
ics of atoms in a symmetric double well potential, where 
the atoms are in an initially fragmented state was also 
analyzed in (Tt| . 

Bose Einstein condensates in triple well structures have 
been analyzed and the stationary solutions in the mean 
field approximation were found in Ref. |l8[ . Three well 
systems show chaotic solutions and the dynamics of 
atoms in a three well potential is sensitive to the initial 
conditions of the system This means that one can 
control the dynamics of the system not only by vary- 
ing the physical parameters of the potential, but also by 
changing the initial conditions. 

The authors of Ref. |2fJ have recently proposed a single 
atom transistor in a ID optical lattice. A quantum inter- 
ference phenomenon is used to switch the flux of atoms 
in a lattice through a cite containing a single impurity 
atom. Finally, Ref. [2lJ discusses "atomtronic" diodes 
and transistors which are direct analogs of their electronic 
counterparts. These devices use cold atoms in an optical 
lattice instead of electrons in a crystal. With the addition 
of impurities both N type and P type "semiconductors" 





FIG. 1: The geometry of a BEC transistor. When the num- 
ber of atoms in the middle well is small, tunneling from the 
left into the right well is negligible (a). This is due to the fact 
that the chemical potential of the middle well does not match 
that of the two other wells (c). Placing atoms in the middle 
well increases the chemical potential due to interatomic inter- 
actions (d) and enables tunneling then atoms tunnel from the 
left into the right well. This happpens because atom-atom 
interactions increase the energy of the middle guide (b). 



may be constructed. A bipolar junction transistor could 
be then built with a NPN or PNP sandwich. 

The rest of the paper is organized as follows. SectioiilTTl 
contains derivation of the general equations of motion for 
a BEC in an n-well potential with arbitrary shapes of the 
wells and discussion of the limits of validity of the model. 
In Sec. I1I1I we specialize our discussion to the case of a 
three well structure. Sec. lIII A"l is devoted to the analysis 
of the equations of motion in the mean field limit and in 
Sec. IIII Bl we will compare the results of the mean field 
to a second quantization calculation. Section lTVI contains 
estimates of the physical parameters for the device and 
discusses the possibility of its experimental realization. 



II. EQUATIONS OF MOTION 

The Hamiltonian for a system of interacting bosons in 
an external potential V(x) is of the form 



H = / dx¥ 



-— v 2 + v> 

Im 



M> + ^ / dx¥¥^. 

(1) 

Here ^! is the field operator and U§ — iira s h 2 /m, where 
m is the atomic mass and a s is the s-wave scattering 



length. For notational simplicity we are considering one- 
dimensional case. Extension to two or three dimensions 
is straightforward. 

In the standard basis of eigenfunctions ipi of the linear 
part of the Hamiltonian 



2m 



V(x) 



ipi = Mliipi 



the field operator is represented as 



(2) 



(3) 



where at is the destruction operator for the mode ipi. 
These operators satisfy the canonical commutation rela- 
tions 



[a t ,a]} = Sij, 
[a;, aj] = 0. 



(4) 



The potential V(x) consists of n weakly coupled potential 
wells. The eigenmodes ipi are "nonlocal" and extend over 
several potential wells. As discussed above, we are inter- 
ested in calculating the number of atoms in each well as 
a function of time. A more convenient basis in this case 
corresponds to a set of modes (pi localized in each poten- 
tial well with the corresponding destruction operators bi 
so that 



(5) 



The operators b are linear superpositions of the operators 
a, 



E 



(6) 



where u is the transformation matrix determined by the 
condition of localization of the modes (pi. 

Requiring that the destruction operators bi satisfy the 
canonical commutation relations identical to those of 
Eq. (D 



[bi,b]} = 

Mi] = o, 



(7) 



implies the unitarity of the transformation matrix u: 
J2 m u mi u mj = F° r bound states all modes ipi can 
be chosen real and the transformation matrix u can be 
chosen real and orthogonal. 

The transformation from the " nonlocal" basis ipi to the 
"local" basis <pi is given by the relations 




FIG. 2: A schematic of a multi-well non-symmetric potential 
structure with two adjacent wells shown. The points Xk-i, x k 
and Xk+i are chosen between the wells where the eigenmodes 
ipk are exponentially small. 



The operators bi are associated with the local modes of 
the n-well structure. For the purposes of the subsequent 
analysis we will need to know only the lowest local mode 
in each potential well. It means that there are n local 
modes <pi and the coefficients tiy should be chosen so 
that the function (pi be localized in the i-th potential 
well. 

To quantify the degree of localization, we set points 
Xo,Xi, ...,X n somewhere between the wells where ampli- 
tudes of the modes ipk are exponentially small. This pro- 
cedure is shown schematically in Fig. [21 The degree of 
localization of the mode <pk in the A;-th well is character- 
ized by the localization parameter 



fk 



dx\(f>k 



(9) 



The localization parameter fk can be rewritten in terms 
of the transformation matrix u as 



fk = ^2u mk Ui k J k (l,m), 

III .1 



(10) 



where the matrix Jk(l,m) is given by the expression 



Jk(l,n%) = / dxip* m ipi. 

J Xk-l 

To localize the modes d>k, we maximize the function 



(11) 



3 



(8) 



/ = £/fc = ^2,u mk u* lk J k (l,m) (12) 

k krnl 



4 



subject to constraints 

^ ] UmkUif. = S. n 



(13) 



The maximization results in the set of n 2 equations 

J~] [Jj(i, m) - Xj m ] u mj = 0, (14) 

m 

where 1 < i, j < n and where Xi m — A^ are Lagrangian 
multipliers. 

The set of equations (|14f) can be written in a more 
transparent form as a set of n matrix equations (j = 
1,2,. ..n) 



(jj - AJ \uj) 



0. 



(15) 



where Jj and A are Hcrmitian matrices with the elements 

i is the col- 
1,2,. ..n). 



Jj(i,m) and Ai m , respectively, and where |'i 
umn vector of u with the elements u m j (m 
The equation of constraints, Eq. (| 1 3|) . becomes 

(Ui\uj) = Sij. 



(16) 



In the limit of negligibly small coupling between the wells, 
the column vectors \uk) of the transformation matrix u 
are exact eigenvectors of the operators Jj because the 
latter in this limit reduce to Jj — \uj)(uj\. The matrix 
of Lagrange multipliers in this limit becomes the iden- 
tity matrix. This observation suggests that for nonzero 
coupling between the wells the vectors \uk) can be found 
perturbatively starting from the eigenvectors \wk) of Jj, 
with eigenvalues close to one: 



Jk\w k ) = Hk\w k ) = (1 - ea k )\w k ), 



(17) 



where e <C 1 characterizes relative coupling strength be- 
tween the wells. The eigenvectors \wk) form a nonorthog- 
onal basis set with (wi\wj) — 0(e) for i ^ j. 

Solution of Eq. (|17|l in the \uk) basis to the first order 
in e yields 

K) = \u k ) + J2\ u o)( u j\Jk\uk), (18) 

where we have used the fact that Jk — \uk)(uk\ = O(e). 
Inversion of Eq. H18|l yields 



\ui) = \Wi) -^2\Wj){uj\Ji\ 



(19) 



Using the orthogonality conditions for \ut} up to the first 
order in e and the condition (uj| Jj\uj) — (ui\ Ji\Uj) that 
follows from Eq. (|15|L results in the relation 



(Ui\Jj\Uj) = -{w^w-j) 



(20) 
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FIG. 3: An example of calculation of local modes. Graphs 
(a)-(c) show the three lowest global eigenmodes of the po- 
tential V (dotted line). These eigenmodes are nonlocal with 
large probability density in two or more wells. Graphs (d)- 
(f) show local modes, which are linear combinations of the 
nonlocal eigenmodes. 



yielding the final expression for the vectors \ui) in terms 
of \wi): 



\Ui) = K) - ~^2\Wj)(Wj\Wi). 

1 i*i 



(21) 



To calculate the local modes, one thus finds eigenvectors 
\wi) of Ji with eigenvalues close to one for i = 1, 2, • • • n. 
The columns of the transformation matrix are then given 
by Eq. I|21|) . The local modes are found using Eq. JSJ. 
An example of such calculation is shown in Fig. |3J 

Overall signs of the local modes 4>i in Eq. ijSJ are ar- 
bitrary being determined by sign choices for the global 
modes ipi . These signs can be changed if needed because 
if 4>i is a local eigenmode, so is — 4>i. Changing the sign 
of (f>i amounts to changing the sign of the z-th row of 
the transformation matrix u which leaves it unitary. To 
be able to unambiguously specify the value of the rela- 
tive phase shift between condensates in different potential 



5 



wells, the overall signs of the local modes <pi will be fixed 
by requiring that each cigcnmodc 4>i be positive in the 
region of its localization between Xi and a;»+i. 

In terms of the destruction operators bi of the local 
modes the Hamiltonian Eq. (JTJ can be written as 

H = J2 MlkUkiut/fy + y E ^ ( b lY b l ( 22 ) 

where fii is the eigenfrcqucncy of the i-th mode ipi given 
by Eq. J2J), Uij is the transformation matrix and Ki is the 
overlap integral 

m = J dx\&\ 4 . (23) 

The equations of motion for the operators bi in the 
Heisenberg picture are given by 

ifi-rbi = t&luuuutfy + J7 (24) 

jk 

The diagonal terms 

k l =Y,^k\u kl \ 2 (25) 

k 

111 Eq. Jill) have the meaning of eigenfrequencies of the 
local eigenmodes in the absence of coupling between the 
wells and the nondiagonal terms 

Aij = ^ tt k u ki u* kl (26) 

are coupling frequencies between the i-th and j-th wells. 
Since for bound states u can be chosen real, the matrix 
of the coupling frequencies is real and symmetric: Ay = 
Aji. The coupling strength is exponentially dependent 
on the distance between the wells and usually only the 
nearest- neighbor coupling should be taken into account. 

III. THREE- WELL STRUCTURE 

In the following we shall specialize our analysis to the 
case of a potential consisting of three potential wells. 
These will be referred to as the left, middle and right 
well, respectively. The left well serves as a source of 
atoms. The number of atoms Ni in this well is kept 
nearly constant and is considerably larger than the num- 
ber of atoms initially placed and subsequently tunneling 
into the middle or the right wells. The dynamics in the 
left well is therefore unaffected by that in the other two 
wells. This dynamics is factored out and the destruction 
operator for the left well bi is replaced by a c-number: 
bi =-> v 7 ^- 

The Hamiltonian Eq. (|22ll reduces to 

H = h{A m - ^b^bra + h{A r - ^)b\b r 
+H&lmvNlb m + A mr b^ m b r + h.c.) 

+ Y^(bl) 2 b 2 m + ^n r (bi) 2 bl (27) 



where Aj, and Ki are given by Eqs. (|2*5|l . i(2!)|) and 
J2HJ, respectively, h.c. means Hermitian conjugate and 
fi = hA l + kiUqNi. 

As discussed at the end of Sec. |D] the overall sign of 
the local modes & has been fixed by requiring that they 
be positive in the region of their localization. With this 
choice, the coupling frequencies Ay between different 
wells (see Eq. (J2SJ0 are negative. This is easily ascer- 
tained using the simplest example of a symmetric two- 
well structure where the two local modes are proportional 
to a sum and a difference of the two global modes. Nor- 
malizing the Hamiltonian Eq. (|27() to the positive energy 
— hA mr — h\A mr \ brings it to its final dimensionless form 

— — 1 = 0;^^^ + LJ r blb r - (Db m + b^br + h.C.) 

+ %^U 2 ^ + f(4) 2 ^ (28) 

where Z t = -U Ki/hA mr , u)i = {n — Aj)/A mr and 
D = A lm vHVi/A mr . 

The Heisenberg equations of motion for the destruction 
operators bi and b r in the dimensionless variables 
take the form 

i-^-frm = (w TO + Z m b\ n b m )b m - D-b r 
i^-b r = (uj r + Z r b\b r )b r — b m , (29) 

(XT 

where the dimensionless time r is given by the relation 
t | A mr \t. 

A. Mean-field 

In this section we shall present results of analysis of 
Eq. (EH! in the mean-field limit corresponding to rela- 
tively large atomic populations in all wells, when the op- 
erators b m and b r can be treated as complex numbers. 

Figure ^demonstrates control of atomic population in 
the right well by population in the middle well with the 
absolute gain that is considerably larger than one. Pa- 
rameters for Fig.^are ui m = —1.3, uj r = 0.5, Z m D 2 = 1 
and Z r D 2 = 0. The right well is initially empty, b r (0) = 
0. Parameters of the wells are chosen so that if no atoms 
are initially placed in the middle well (b m (Q) — 0), the 
tunneling from the source (the left well) to the middle 
well is strongly suppressed and the population in the right 
well remains at a low level. This situation is illustrated 
by a dotted line in Fig. 

Placing some number of atoms in the middle well re- 
sults in a much larger tunneling rate from the left to the 
right well through the middle well as shown by a solid 
line corresponding to the initial condition 6 m (0) = D. 
The increase in the tunneling rate can be observed for a 
range of values of the relative phase of the condensates 
in the left and middle wells. The dashed curve obtained 
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FIG. 4: The number of atoms in the right well as a function 
of interaction time for different initial number of atoms in the 
middle well. The dotted curve corresponds to initially empty 
middle well, 6 m (0) = 0, and the solid curve to 6 m (0) = D. 
The dashed curve corresponds to the initial condition b m (0) = 
Dexp(in/2). For all curves Lo m = —1.3, ui r — 0.5, Z m D 2 — 1 
and Z r D 2 = 0. 



for the initial condition 6 m (0) = -Dexp(i7r/2), i.e., cor- 
responding to the 7r/2 relative phase shift between the 
condensates in the left and middle wells, exhibits quali- 
tatively similar behavior. Note that the output number 
of atoms in the right well (r = 20) is about 50 — 60 times 
larger than the input number of atoms in the middle well. 
In other words, the output number of atoms in the right 
is controlled by that in the middle well with the absolute 
gain G = N rt0ut /Ni tin ra 50 - 60. 

Populations in the middle and right wells as functions 
of the interaction time are shown if Fig.|S]for b m (0) = D. 
All other parameters are the same as in previous graphs. 
The solid curve is the population of the right well and 
the dashed curve is the population of the middle well. 
Figure demonstrates that the population of the middle 
well stays about an order of magnitude below that for 
the right well. The middle well serves a gate controlling 
the rate of atomic flow from the source to the right well. 
The atoms tunneling from the source to the right well 
pass through the middle well without being accumulated 
there. 

The output number of atoms in the right well as a 
function of the input number of atoms in the middle 
well is shown in Fig. Parameters for this figure 
are the same as for Fig. 0J i.e., Lu m = —1.3, uj r = 0.5, 
Z m D 2 = 1 and Z r D 2 = 0. The solid curve corresponds 
to the zero initial phase shift between the condensates 
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FIG. 5: The number of atoms in the middle (dashed) and 
right (solid curve) well as a function of interaction time for 
Lo m = -1.3, Lj r = 0.5, Z m D 2 = 1, Z r D 2 = and b m = D 



in the left and middle wells and the dotted curve to 
the 7r/2 shift. This figure demonstrates rapid switching 
from small to large tunneling rates in the region around 
\b m ,in\ 2 D 2 ~ 0.5 with subsequent saturation at the level 
G = N r , ou t/N m ,in « 50 — 60. In the switching region, a 
small change in the population of the atoms in the mid- 
dle well results in a large difference in the population in 
the right well. 

Figure |7| shows output population in the right well 
(r = 20) as a function of the modal frequency oj m of 
the middle guide for different values of the input number 
of atoms in the middle well. This figure demonstrates 
switching for different values of the number of atoms ini- 
tially in the middle well. The dotted line corresponds 
to initially empty middle well. For this curve, the max- 
imum tunneling rate corresponds to the region around 
ui m = —0.5. If the frequency of the middle well is low- 
ered beyond this value, the number of atoms that tunnel 
into the right well becomes small. The solid curve corre- 
sponds to the initial condition b m (0) = D. This curve is 
qualitatively similar to that for an initially empty mid- 
dle well, but the maximum has moved to a lower value of 
uj m = —1.3. The dashed curve corresponds to the initial 
condition b m (0) — D exp(i7r/2). This curve still has a 
maximum around u) m = —1.3. For this value of u> m , the 
number of atoms that tunnel into the right well when 
the initial atoms have either zero or 7r/2 phase shift is 
about the same. For the initially empty middle well, the 
population in the right well remains small. 

As opposed to an electronic transistor, the amplifica- 
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FIG. 6: The output number of atoms in the right well (r = 
20) as a function of the number of atoms initially placed in 
the middle well. 



80 




FIG. 7: The output number of atoms in the right well (r = 
20) as a function of the frequency of the middle well iu m - The 
dotted line corresponds to the initially empty middle well, 
bm (0) = 0, the solid line to b m = D, and the dashed line to 
b m = -Dexp(i7r/2). 




FIG. 8: The output number of atoms in the right well as 
a function of the relative phase of the atoms placed in the 
middle well. 

tion and switching in the three- well structure is a coher- 
ent effect and depends on the relative phase between the 
condensates in the left and middle wells. To investigate 
sensitivity of the previously obtained results to the value 
of the relative phase angle, we kept the input number 
of atoms in the middle well fixed at |^4 m (0)/£>| 2 = 1.2 
and changed the relative phase angle. The results are 
given by Fig. |H1 showing the output number of atoms in 
the right well (r = 20) as a function of the phase angle. 
All other parameters are the same as in previous figures. 
Figure |S] demonstrates that in the amplification regime 
the number of atoms that tunnel into the right well is 
nearly independent of the initial phase angle, as long as 
this angle is roughly in the range between —tt/2 and 7r/2. 

B. Second-quantization results 

This section presents results of analysis of the 
three-well structure in the framework of the second- 
quantization formalism. This allows us to estimate the 
region of applicability of the mean-field approach of sec- 
tion llll XI evaluate intrinsic quantum-mechanical uncer- 
tainty due to finite number of atoms and extend previous 
results to the limit of small number of atoms. 

In the dimensionless variables, the state vector of the 
system \ip(t)) evolves according to the equation 

i±\^)=H\^ (30) 
where H is the second-quantized Hamiltonian given by 



8 



Eq. (J2BJ. The state vector can be represented in terms 
of the joint number states \n m , n r ) as 



>.:) 



(31) 



with the decomposition coefficients given by c m . r = 
(n m ,n r \ip)- Equation is transformed to the set of 
ordinary differential equations that describe the evolu- 
tion of the decomposition coefficients 



i,j,k,l 



(32) 



In simulations, the set of Eqs. ll-'il'l) has been truncated 
by keeping only the values of n m and n r such that n r + 
Ti m < N max . The value ofN max was chosen so that N max 
was several times larger than the sum (n r ) + (n m ). 

Initial conditions for the system of equations Eq. i|32|) 
corresponded to zero initial number of atoms in the right 
well with the atoms in the middle well being in a coherent 
state: 



M0)> =e"W ^M>- 

71 — V 



(33) 



Here the complex parameter a is given by a = 
V (-/V m )(0)e* v , where (iV m )(0) is the average number of 
atoms initially placed in the middle well and tp is the 
phase difference between the atoms in the middle and 
left wells. 

The transition to the mean-field limit corresponds to 
increasing the input number of atoms in the middle well 
(N m )(0) while keeping the ratio ui m /Z m (N m }(0) con- 
stant. Equations 1|28[) show that the results of action 
of the destruction operators on the state vector scales 
as D provided the parameters Z m D 2 and Z r D 2 are kept 
constant. Thus, the transition to the mean-field limit can 
be implemented by setting the initial number of atoms 
in the middle well proportional to D 2 and increasing the 
value of coupling D between the left and middle wells 
while keeping the parameters Z m D 2 and Z r D 2 constant. 

The average number of atoms in the right well (N r ) as 
a function of the interaction time is shown in Fig. |^| for 
three different values of D 2 . The parameters for this fig- 
ure are u> r = 1, u rn = —0.5, Z m D 2 — 1/4, Z r D 2 = 
and \a\ 2 = (N m )(0) = D 2 . The phase angle of the 
coherent state is zero. The solid line is the mean-field 
limit. The dotted, dash-dotted and dashed lines corre- 
spond to D 2 = I, D 2 = 4 and D 2 = 8, respectively. 
Figure demonstrates good convergence of the second- 
quantization results to the mean field limit as D 2 is in- 
creased. The D 2 — 1 curve deviates from the mean-field 
limit for large value of r, but all other curves lie progres- 
sively closer to the mean-field curve as the parameter D 
increases. 

The output (t = 10) probability distribution P(N r ) of 
finding N r atoms in the right well versus N r is shown if 
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FIG. 9: The average number of atoms in the right well 
as a function of the interaction time for three different val- 
ues of the coupling between the left and the middle wells. 
Atoms in the middle well are initially in a coherent state with 
(Nm}(0) = D 2 and zero phase. The solid line is the result of 
the mean-field calculation, the dotted line is the result of the 
second-quantization calculation with D 2 = 1, the dashed line 
corresponds to D 2 = 4, and the dash-dotted one to D 2 — 8. 




5 10 15 20 25 30 35 40 45 50 
N/D^ 



FIG. 10: The probability P(N r ) of finding N r atoms in the 
right well (r = 10). Atoms in the middle well are initially in a 
coherent state with a — D. The dash-dotted line corresponds 
to D 2 — 1, the dashed line to D 2 = 4 and the solid line to 
D 2 = 8. 
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Fig. The dash-dotted line corresponds to D 2 = 1, 
the dashed line to D 2 — 4 and the sold line to D 2 = 8. 
Since the number of output atoms scales as D 2 , the hor- 
izontal axes is scaled as N r /D 2 to keep position of the 
maximum and the width of the curves more of less the 
same for different values of D 2 . As a result, the ver- 
tical axes shows not P(N r ), but the product P(N r )D 2 
to keep the height of the curves approximately the same 
for different values of D 2 . The total "area under the 
curve" (strictly speaking it is a sum, not an integral) for 
all curves is equal to one. The dash-dotted curve cor- 
responding to D 2 = 1 shows bimodal distribution with 
a relatively large probability of finding atoms near zero 
N r in addition to the main peak near N r /D 2 i=s 31, the 
latter being very close to the mean-field result. The dif- 
ference between the mean-field and second-quantization 
results, previously seen in Fig.Elfor D 2 = 1, is due to the 
part of the probability distribution near zero that pulls 
down the average. As the coupling D to the source is 
increased, only the single-humped part of the probability 
P(N r ) centered at the mean-field result remains. The 
output relative standard deviation AN r / (N r ) is equal to 
0.35, 0.08 and 0.04 for D 2 = 1,4 and 8, respectively. 

Comparison of the mean-field and second-quantization 
results carried out for the same parameters as above but 
the relative phase angle between the condensates equal 
to (f> = 7r/2 yielded conclusions very similar to those sum- 
marized by Figs. 151 and 1101 

Figures ITTlandlT^lparallel analysis of Figs . I§1 and 1 1 01 for 
the case when the middle well is initially empty, | ck| 2 = 
(N m }(0) — 0. These figures are aimed at verifying that 
the rapid switching from low to high-amplification regime 
predicted by the theory in the mean-field limit can be also 
realized with only few controlling atoms. 

Figure 1111 demonstrates convergence of the second- 
quantization results to the mean-field limit for (N m ) (0) = 
as D 2 is increased. This convergence is similar to 
that shown in Fig. [5] except in this case the second- 
quantization approach gives values somewhat larger than 
the mean- field limit. The reason is explained by Fig. 1121 
which shows the probability P(N r ) of finding N r atoms 
in the right well at r = 10. The probability P(N r ) has a 
pronounced spike at low values of N r . Another noticeable 
feature of Fig. ^| is a wide, nearly flat pedestal extend- 
ing from small values of N r to a maximum value that is 
about twice larger than the average (cf. Fig. Illfl . The 
maximum value slightly decreases as D increases. This 
explains why the second-quantization results are larger 
than the mean-field results. The one-humped shape of 
P(N r ) in Fig. 1101 means that the uncertainty in the out- 
put number of atoms in the high- amplification regime 
is small for even a few controlling atoms in the middle 
well. Figures ^] an( i 133 show that the low-amplification 
region is characterized by both low average number of 
output atoms and by large uncertainty corresponding to 
the average. Indeed, the output relative standard devia- 
tion AN r /(N r ) for the results of Fig. ^| is equal to 1.3, 
1.0 and 0.9 for D 2 = 1,4 and 8, respectively. These re- 
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FIG. 11: The average number of atoms in the right well as 
a function of the interaction time for three different values 
of the coupling between the left and the middle wells. The 
middle well is initially empty. The solid curve is the result 
of the mean-field calculation, the dotted curve is the result of 
the second-quantization calculation with D 2 — 1, the dashed 
curve corresponds to D 2 = 4, and the dash-dotted curve to 
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FIG. 12: The probability P(N r ) of finding N r atoms in the 
right well (r = 10). The middle well is initially empty. The 
dash-dotted line corresponds to D 2 = 1, the dashed line to 
D 2 = 4 and the solid line to D 2 = 8. 
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suits are in contrast to those for the large-amplification 
regime of Fig. 1101 where the standard deviation rapidly 
goes down as the parameter D 2 increases. 



IV. DISCUSSION 



The analysis of Sec. IIIII demonstrates that a Bose 
Einstein condensate in a three well potential shows 
transistor-like behavior with the middle well acting as 
a gate controlling the flux of atoms from the source to 
the drain. In this section we present estimates of the 
characteristic tunneling time for a trapped atom tran- 
sistor and discuss possible gain in the total number of 
atoms. The analysis will be extended to the case of a 
waveguide device, where estimates will be presented for 
the tunneling time, the length of the device and the gain 
in the output flux of atoms. Finally, we summarize the 
results obtained. 



A. Trapped atom transistor 

The parameter D in Eq. (|29|) characterizes the strength 
of coupling of the source (left well) to the gate (medium 
well). It is reasonable to expect that the operational 
parameters of the BEC transistor are such that the con- 
tributions of the nonlinear and linear terms in Eq. I|29|) 
are of the same order of magnitude, i.e. Z m D 2 ss 1. 

The growth curve, shown if Fig. [S] shows the final 
population of the right well as a function of the initial 
population of the middle well, where the atoms are held 
in the traps for a dimensionless time r = 20. This fig- 
ure demonstrates that a change in the population of the 
middle well from N m = OAD 2 to N m = 0.8D 2 results 
in a change in the final population of the right well from 
N r S3 10D 2 to N r S3 60£> 2 . The maximum number of 
atoms that tunnel into the right occurs when the num- 
ber of atoms initially in the middle well is N m sa D 2 . We 
will refer to this number as the saturation number. For 
example, if we take D 2 = 10, a change from 4 to 8 atoms 
in the middle well results in a change from 100 to 600 
atoms in the right well. 

Assume that the potential energy of the middle well is 
a cigar shaped potential of the form 



V(r±,z) = -m {oJ 2 ± rl + uo 2 z z 2 ) 



(34) 



where r± is the coordinate in the radial direction, and 
z is the coordinate in the axial direction. For this po- 
tential, the overlap integral given by the Eq. (|2*3|) can be 
evaluated as 



1 



1 



(2tt) 3 / 2 a\a z 



(35) 



where a± = y/h/mu>_i and a z — ^Jh/raui z are the har- 
monic oscillator lengths. The nonlinearity parameter Z m 



in Eq. (|28|l is given by the expression 

U()K m 



Z m — 



H\A mr \ 



(36) 



where U = Ana s h 2 /m. Using Eq. O in Eq. iJSBJ 
allows one to express the tunneling frequency between 
the middle and the right wells as 



A mr — 



N„ 



a z 



(37) 



where we have also used Z m N m S3 1 to eliminate Z m in 
favor of N m . 

If the middle well is a spherical trap with lu z = cjj_ = 
2ir x 10 3 Hz and the saturation number is D 2 = 10 the 
tunneling frequency between the middle and right well is 



7r x I0 2 rad/sec. 



(38) 



The dimensional time that it takes for atoms to tunnel 
from the left to the right wells is 



t « 2 x 10 _1 sec. 



(39) 



In other words, for the parameters chosen a trapped 
atom transistor can distinguish between 4 and 8 atoms 
in the gate with the characteristic operational time of 
10 sec. This time can be decreased either by increasing 
the frequency of the trap or increasing the value of the 
saturation number. 



B. Waveguide transistor 

In a waveguide transistor the potential wells of Fig. 2] 
are the three guides that run parallel to each other for the 
distance L. The interaction time T = L/v is determined 
by the speed of flow v of the BEC in the guides. The 
field operator for this configuration can be expressed as 



\&(r, t) — exp(ik p z — iuj p t)if>(r,t), 



(40) 



where k p and to p = hk 2 /2m and the carrier wave number 

and frequency, respectively, and ip is the field-operator 
envelope. 

The Heisenberg equations of motion for the field oper- 
ator tp in the co-propagating frame t' = t, z' = z — vt is 
of the form 



2m 



d 2 



(41) 

where v = hk p /m is the velocity of the condensate and 
the primes have been omitted. 

Changes in density as the condensate propagates 
through the transistor occur at a length scale Lbec- 
We assume that the kinetic energy associated with the 
longitudinal direction is small in comparison with the 
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characteristic energy Ml associated with the transverse 
eigenmodes of the transistor 



Ml > 



2mL 2 BEC 



(42) 



Next, we require that Lbec does not change appreciably 
during the time interval L/v that it takes the condensate 
to propagate through the transistor, 



d~t LBEC 



(43) 



With Eqs. (g2J and fulfilled, the dispersive term 
(d 2 jdz 2 ) in Eq. (|41|) can be neglected and the coordi- 
nate z becomes a parameter. Propagation of different 
"slices" of the condensate (parametrized by the coordi- 
nate z) through the transistor can be analyzed indepen- 
dently. 

Represent the field operator ifj as 



4> = ^2M r ±)h(z,t), 



(44) 



where <fr(r±) is the i-th transverse local mode and bi(z, t) 
is the destruction operator that destroys an atom in the 
i-th local mode at the coordinate z. Note that bi now 
has dimension of m" 1 / 2 . 

Using Eq. (|44|) in Eq. (|41|) with the dispersive term 
dropped, results in the equations of motion in the Heisen- 
berg picture for the operators bi that are of the same form 
as Eq. (EH): 



ih—bi(z,t) = ^2 hQ k UkiU* kj bj + UoKib\bi, 



(45) 



As in the case of a trapped device, the left guide will 
be treated as a reservoir of atoms corresponding to the 
replacement bi — ► y/ni, where n; is the density of atoms 
(number of atoms per unit length). The equations of 
motion for the atoms in the middle and right guide, in 
dimensionless form, become 



*^T 6 ™( Z ' = ( w ™ + z mb\nb m )b m — D — b r 



i—b r (z,t) — (uj r + Z r b\b r )b r 
dr 



(46) 



where the dimensionless parameters are Zi — 
-U Ki/hLA mr , Ui = (m - Ai)/A mn D = 
AlmV Lni I A mr , and L is the length of the transis- 
tor. The destruction operators are normalized to 
K = bi/yL, and the primes have been dropped. Since 
the equations for each "slice" in z are of the same form 
as Eq. (|29[) . the analysis of Sec. IIIII is valid for each 
"slice" separately. 

As with the case of a trapped atom transistor, we take 
Z m D 2 = 1 and use the fact that the largest tunneling 
rate corresponds to n m w D 2 . We refer to this as the 
saturation density, since n m is the normalized density of 



atoms and not the total number as it was with a trapped 
atom device. Next, we assume that the middle waveguide 
can be described by the potential 



V(r ± ) = \mu) 



(47) 



where uj± is the transverse frequency of the guide. The 
overlap integral associated with this potential is 



1 1 

27 



(48) 



where a± is the transverse oscillator length and aj_ = 
\JhJrrm±. Using Eq. (|3*6*)) and l|48|l. we evaluate the 
coupling frequency between the middle and right guides 
as 



I A I ^ " m 

|A mr | ~ a s uj±—j- 



(49) 



In terms of the velocity of the atoms v and the flux en- 
tering the middle guide <!> m the density can be expressed 
as n m /L — § m /v, and Eq. (|49|l takes the form 



|A„ 



* r. 



a s uj±- 



(50) 



Assuming that the guide has a transverse frequency of 
10 kHz, the velocity 5 cm/sec and the saturation flux is 
10 5 atoms/sec, we can evaluate the coupling frequency 
between the right and middle guide as 



A, 



2n x 10 2 rad/sec. 



The dimensional switching time is 

t w 2 x 10 _1 sec, 
and the length of the device is 



lcm. 



(51) 



(52) 



(53) 



This length can be decreased by slowing the velocity of 
the atoms, increasing the saturation flux or increasing 
the transverse frequency of the waveguide. 

With the above numbers, a change in the input flux 
of the middle guide from 0.4 x 10 5 atoms/sec to 0.8 x 
10 5 atoms/sec results in a change of flux in the out- 
put of the right guide from 10 6 atoms/sec to about 
10 7 atoms/sec. 

To summarize, we have presented a theoretical anal- 
ysis of a Bose Einstein condensate in a nonsymmetric 
three- well potential which shows transistor-like behavior. 
We demonstrated the control of atomic population in the 
right well by the population in the middle well with an 
absolute and differential gains considerably larger than 
one. The second-quantization formalism was then used 
to evaluate the quantum-mechanical uncertainty due to a 
finite number of atoms and extend the mean-field results 
to the limit of a small number of atoms. 

The BEC transistor can turn out to be useful in pre- 
cision measurements. The number of atoms that tunnel 
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from the source to the drain is very sensitive to the num- 
ber of atoms in the gate. This fact can be used to de- 
tect and amplify small changes in the number of atoms 
in the gate. A waveguide based transistor is capable of 
operating continuously and can be used to measure time- 
dependent phenomena. Applications of this device may 
include measurement of inertial changes and electromag- 
netic fields. It is possible to envision potentially more in- 
teresting applications by combining several such devices 
so that, e.g., the amplified output of the first transistor 
serves as control for the second. 
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